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Abstract 

Functional relation for commuting quantum transfer matrices of quantum integrable models is 
identified with classical Hirota's bilinear difference equation. This equation is equivalent to the 
completely discretized classical 2D Toda lattice with open boundaries. The standard objects of 
quantum integrable models are identified with elements of classical nonlinear integrable difference 
equation. In particular, elliptic solutions of Hirota's equation give complete set of eigenvalues of the 
quantum transfer matrices. Eigenvalues of Baxter's Q-operator are solutions to the auxiliary linear 
problems for classical Hirota's equation. The elliptic solutions relevant to Bethe ansatz are studied. 
The nested Bethe ansatz equations for Afe_i-type models appear as discrete time equations of motions 
for zeros of classical r-functions and Baker- Akhiezer functions. Determinant representations of the 
general solution to bilinear discrete Hirota's equation and a new determinant formula for eigenvalues 
of the quantum transfer matrices are obtained. 
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1 Introduction 



In spite of the diversity of solvable models of quantum field theory and the vast variety of methods, the 
final results display dramatic unification: the spectrum of an integrable theory with a local interaction 
is given by a sum of elementary energies 

E = Y^s{u{), (1.1) 

i 

where Ui obey a system of algebraic or transcendental equations known as Bethe equations ||l[], The 
major ingredients of Bethe equations are determined by the algebraic structure of the problem. A typical 
example of a system of Bethe equations (related to Ai-type models with elliptic i?-matrix) is 
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where a{x) is the Weierstrass a-function and 
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Entries of these equations which encode information of the model are the function e{u) (entering through 
4>{u)), quasiperiods wi, uj2 of the cr-function, parameters rj, v, and size of the system N. 

Different solutions of the Bethe equations correspond to different quantum states of the model. 

In this paper we show that these equations, which are usually considered as a tool inherent to the 
quantum integrability, arise naturally as a result of the solution of entirely classical non-linear discrete 
time integrable equations. This suggests an intriguing interrelation (if not equivalence) between integrable 
quantum field theories and classical soliton equations in discrete time. In forthcoming papers we will 
show that the Bethe equations themselves may be considered as a discrete integrable dynamical system. 

In 1981 Hirota |^ proposed a difference equation which unifies the majority of known continuous 
soliton equations, including their hierarchies Q. 

A particular case of the Hirota equation is a bilinear difference equation for a function r(n, /, m) of 
three discrete variables: 

aT{n,l + l,m)T{n,l,m + l) + l3T{n,l,m)T{n,l + l,m + l) +jT{n+l,l + l,m)T{n — l,l,m + l) — 0, (1.4) 

where it is assumed that a + /3 + 7 = 0. Different continuum limits at different boundary conditions 
then reproduce continuous soliton equations (KP, Toda lattice, etc). On the other hand, T{n,l,m) can 
be identified Q] with the r-function of a continuous hierarchy expressed through special independent 
variables. 

The same equation (with a particular boundary condition) has quite unexpectedly appeared in the 
theory of quantum integrable systems as a fusion relation for the transfer matrix (trace of the quantum 
monodromy matrix). 

The transfer matrix is one of the key objects in the theory of quantum integrable systems |^ . Transfer 
matrices form a commutative family of operators acting in the Hilbert space of a quantum problem. Let 
Ri.Aiu) be the i?-matrix acting in the tensor product of Hilbert spaces Vi (E) K4. Then the transfer matrix 
is a trace over the auxiliary space of the monodromy matrix. The latter being, the matrix product 
of N i?-matrices with a common auxiliary space: 

TA{u\yi) = RN,A{u-yN)---R2.,A{u-y2)Ri,A{u-yi), 

Ta{u) = tr^TU(u|yO . (1.5) 

The transfer matrices commute for all values of the spectral parameter u and different auxiliary spaces: 

[Ta{u),Ta'{u )]={). (1.6) 



2 



They can be diagonalized simultaneously. The family of eigenvalues of the transfer matrix is an object of 
primary interest in an integrable system, since the spectrum of the quantum problem can be expressed 
in terms of eigenvalues of the transfer matrix. 

The transfer matrix corresponding to a given representation in the auxiliary space can be constructed 
out of transfer matrices for some elementary space by means of the fusion procedure 0, Q. The 
fusion procedure is based on the fact that at certain values of the spectral parameter u the i?-matrix 
becomes essentially a projector onto an irreducible representation space. The fusion rules are especially 
simple in the Ai-case. For example, the i?i^i(M)-matrix for two spin-1/2 representations in a certain 
normalization of the spectral parameter is proportional to the projector onto the singlet (spin-0 state) 
at u = +2 and onto the triplet (spin-1 subspace) at u — —2, in accordance with the decomposition 
[1/2] + [1/2] = [0] + [1]. Then the transfer matrix T^{u) with spin-1 auxiliary space is obtained from the 
product of two spin-1/2 monodromy matrices Tl{u) with arguments shifted by 2: 

r2i(M) = tr[i](i?ia(-2)fi(^i + l)Ti(^.-l)i?i,i(-2)). 

A combination of the fusion procedure and the Yang-Baxter equation results in numerous functional 
relations (fusion rules) for the transfer matrix [^, They were recently combined into a universal 
bilinear form [ p^ , [ pT] . The bilinear functional relations have the most simple closed form for the 
models of the ^fc_i-series and representations corresponding to rectangular Yonng diagrams. 

Let T^(u) be the transfer matrix for the rectangular Young diagram of length a and height s. If ry 
can not be represented in the form rj — ritoi + r2UJ2 with rational ri, r2 (below we always assume that 
this is the case; for models with trigonometric i?-matrices this means that the quantum deformation 
parameter q would not be a root of unity), they obey the following bilinear functional relation: 

T^{u + l)T^{u - 1) - T-+MT^-M - T-+\u)Tr\u) . (1.7) 



Since Tgiu) commute at different u, a, s,, the same equation holds for eigenvalues of the transfer 
matrices, so we can (and will) treat T"(u) in eq. (1.7) as number-valued functions. The bilinear fusion 
relations for models related to other Dynkin graphs were suggested in ref. [O]. 



Remarkably, the bilinear fusion relations (1/7) appear to be identical to the Hirota equation (1.4) 
Indeed, one can eliminate the constants a, /?, 7 by the transformation 



r(n, /, m) = 



(1-^7/0;)'" 



T„(/, m). 



so that 



^{l + l,m)T„{l,m+ 1) - Tn{l,m)T„{l + l,m + 1) + t„+i(? + 1, r7i)T„_i(/, m + 1) = 



and then change variables from light- cone coordinates n,l,m to the "direct" variables 



n. 



s = / + m, u — I — m ~ n, 



i{l,m) 



(1.8) 



(1.9) 



At least at a formal level, this transformation provides the equivalence between eqs. ( |1.7| ), ( |1.4| ) and 



In what follows we call eq. ( |1.§| ) (or (1.7)) Hirota's bilinear difference equation (HBDE). 

Leaving aside more fundamental aspects of this "coincidence", we exploit, as a first step, some 
technical advantages it offers. Specifically, we treat the functional relation (1.7) not as an identity but as 
a fundamental equation which (together with particular boundary and analytical conditions) completely 
determines all the eigenvalues of the transfer matrix. The solution to HBDE then appears in the form 
of the Bethe equations. We anticipate that this approach makes it possible to use some specific tools of 
classical integrability and, in particular, the finite gap integration technique. 



The origin of T°(w) as an eigenvalue of the transfer matrix (1.^ ) imposes specific boundary conditions 
and, what is perhaps even more important, requires certain analytic properties of the solutions. As a 
general consequence of the Yang-Baxter equation, the transfer matrices may always be normalized to 
be elliptic polynomials in the spectral parameter, i.e. finite products of Wcierstrass cr-functions (as in 
(1.3)). The problem therefore is stated as of finding elliptic solutions of HBDE. 
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A similar problem appeared in the theory of continuous soliton equations since the works , , 
wherein a remarkable connection between motion of poles of the elliptic solutions to the KdV equation 
and the Calogero-Moser dynamical system was revealed. Elliptic solutions to Kadomtsev-Petviashvili 
(KP), matrix KP equations and the matrix 2D Toda lattice (2DTL) were analyzed in the Refs. [Q, 
| p5[ , [ p^ , respectively. It was shown, in particular, that poles of elliptic solutions to the abelian 2DTL 
(i.e. zeros of corresponding r-functions and Baker- Akhiezer functions) move according to the equations 
of motion for the Ruijsenaars-Schneider (RS) system of particles p^ . 

Analytic properties of solutions to HBDE relevant to the Bethe ansatz suggest a similar interpretation 
of Bethe ansatz equations. We will show that the nested Bethe ansatz for Ak-i-type models is equivalent 
to a chain of Backlund transformations of HBDE. The nested Bethe ansatz equations arise as equations 
of motion for zeros of the Baker- Akhiezer functions in discrete time (discrete time RS system |^ . The 
discrete time variable is identified with level of the nested Bethe ansatz. 

The paper is organized as follows. In Sect. 2 we review general properties and boundary conditions 
of solutions to HBDE that yield eigenvalues of quantum transfer matrices. In Sect. 3 the zero curvature 
representation of HBDE and the auxiliary linear problems are presented. We also discuss the duality re- 
lation between "wave functions" and "potentials" and define Backlund flows on the set of wave functions. 
These flows are important ingredients of the nested Bethe ansatz scheme. For illustrative purposes, in 
Sect. 4, we give a self-contained treatment of the Ai-case, where major part of the construction contains 
familiar objects from the usual Bethe ansatz. Sect. 5 is devoted to the general Ak-i-case. We give a 
general solution to HBDE with the required boundary conditions. This leads to a new type of determi- 
nant formulas for eigenvalues of quantum transfer matrices. The explicit form of generalized Baxter's 
relations (difference equations for Qt{u)) is presented. They are used for examining the equivalence to 
the standard Bethe ansatz results. In Sect. 6 a part of the general theory of elliptic solutions to HBDE 
is given. Sect. 7 contains a discussion of the results. 



2 General properties of solutions to Hirota's equation relevant 
to Bethe ansatz 

2.1 Boundary conditions and analytic properties 



HBDE has many different solutions. Not all of them give eigenvalues of the transfer matrix (1.5). There 



are certain boundary and analytic conditions imposed on the transfer matrix (1.5). 

(i) It is known that T^{u), the transfer matrix in the most antisymmetrical representation in the 
auxiliary space, is a scalar, i.e. it has only one eigenvalue (sometimes called quantum determinant 
detqTs{u) of the monodromy matrix). It depends on the representation in the quantum space of the 
model and is known explicitly. In the simplest case of the vector representation (one-box Young diagram) 
in the quantum space it is p9| : 

fc-ls-l fe-l 

T^{u) = (l){u ~ s ~ k)Y[Yl(l){u + s + k - 21 - 2p - 2)Y[cl){u + s + k ~ 21), (2.1) 

1=0 p=l 1=1 

T^u) = 1. (2.2) 

These values of T^{u) and T^{u) should be considered as boundary conditions. Let us note that they 
obey the discrete Laplace equation: 

T^{u + \)T^{u - 1) = T^^,{u)Tli{u). (2.3) 

This leads to the boundary condition (b.c.) 

T^iu) = as a < and a > /fc (2.4) 

^It should be noted that equations of motion for the discrete time RS system were already written down in the paper 
|l8[ . However, the relation to elliptic solutions of discrete soliton equations and their nested Bethe ansatz interpretation 
were not discussed there. 
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(with this b.c. eq. ( |1.8| ) is known as the discrete two-dimensional Toda molecule equation |20[, an 
integrable discretization of the conformal Toda field theory pfl). 

(ii) The second important condition (which follows, eventually, from the Yang-Baxter equation) is 
that T^(u) has to be an elliptic polynomial in the spectral parameter u. By elliptic polynomial we mean 
essentially a finite product of Weierstrass ct- functions. For models with rational i?-matrix it degenerates 
to a usual polynomial in u. 



To give a more precise formulation of this property, let us note that eq. (1.7) has the gauge invariance 
under a transformation parametrized by four arbitrary functions Xi of variable: 

Ts{u) Xi{a + U + s)x2{a -u + s)x3{a + u- s)x4{a - u - s)T^{u) (2.5) 

These transformations can remove all zeros from the characteristics a ± s ± u — const. We require that 
the remaining part of all T^{u) should be an elliptic (trigonometric, rational) polynomial of one and the 
same degree N, where N is the number of sites on the lattice (see (1.3)). 

One can formulate this condition in a gauge invariant form by introducing the gauge invariant com- 
bination 



n+\u)Tr\u) 



y:i-) - ■ (2.6) 



We require Y^{u) to be an elliptic function having 2N zeros and 2N poles in the fundamental domain. 
This implies that T^{u) has the general formQ 

N 

T:{u) = Ale^^'^^^^ n '^('K" - ^f''^)) , (2.7) 
where zj"'"-*, fj,(a,s) do not depend on u and the following constraints hold: 

N N 
^^^(a,.+ l) ^ ^(a,s-l)^ ^ ^ ^ ^2.8) 

^(a, s + 1) + /i(a, s - 1) = ^{a + l,s) + fi{a - 1, s) . (2.9) 
Another gauge invariant combination, 

T-{u + l)T-{u-l ) 

is also convenient. 

As a reference, we point out gauge invariant forms of HBDE ||ll| 

(i + n°+i(u))(i + y°_i(u)) 
(i + (y/+i(^.))-i)(i + (rr^(u))-i) 



= -^A_Z_^^^V2 ^ _i _ y:{u) , (2.10) 



Y;(u + i)Y;iu-i) = ^ , , (2.11) 



It can be shown that the minimal polynomial appears in the gauge 

Ca-ls-l a-i \ ^-"^ 

Y\Y\(j3{u + s + a-2l-2p-2)W(j){u + s + a-2l)\ , (2.13) 
;=o p=i 1=1 ) 

where all the "trivial" zeros (common for all the eigenvalues) of the transfer matrix are removed (see 
e.g. 1^). The boundary values at a = 0, /c then become: 

TO(u) = 0(^-1-5), 

T^(u) = 4>{u-s-k). (2.14) 



^This differs from a more traditional expression in terms of Jacobi 0-functions by a simple normalization factor. 
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^From now on we adopt this normalization. 



(iii) The analyticity conditions and b.c. (2.14) lead to a particular "initial condition" in s. It is 
convenient, however, to take advantage of it before the actual derivation. The condition reads 



T^{u) = for any 
This is consistent with (pT7|), (2.14) and implies 



k < s <0, < a < k. 



(u — a) 



(2.15) 



(2.16) 



for < a < fc. 



Under the anal ytici ty conditions (i) and the b.c. (2.14) (and their consequences (2.15), (2.16)) each 
solution to HBDE (1/7) corresponds to an eigenstate of the Afe_i-transfer matrix. 

The same conditions are valid for higher representations of the quantum space. However, in that 
case there are certain constraints on zeros of 4>{u) (they should form "strings"), whence Tf(u) acquires 
extra "trivial" zeros. Here we do not address this question. 



2.2 Pliicker relations and determinant representations of solutions 

Classical integrable equations in Hirota's bilinear form are known to be naturally connected 

p6| , with geometry of Grassmann's manifolds (grassmannians) (see |^ ||2^), in general of an in finite 



dimension. Type of the grassmannian is specified by boundary conditions. Remarkably, the b.c. (2.4) 
required for Bethe ansatz solutions corresponds to finite dimensional grassmannians. This connection 
suggests a simple way to write down a general solution in terms of determinants and to transmit the 
problem to the boundary conditions. Numerous determinant formulas may be obtained in this way. 

The grassmannian Gl^\ is a collection of all (r + l)-dimensional linear subspaces of the complex 
(n + l)-dimensional vector space C""*"^. In particular, G^_|_]^ is the complex projective space P". Let 

X G Gl^i be such a (r+l)-dimensional subspace spanned by vectors x'^' = J2i=o s^P'e', j = 1, . . . , r+1, 
where e' are basis vectors in C"+^. The collection of their coordinates form a rectangular (n+1) x (r+1)- 
matrix Let us consider its (r + 1) x (r + 1) minors 

det(a;|'^) = (?o,«i, • ■ ■ , v), p, g = 0, 1, . . . , r , (2.17) 

pq P 

obtained by choosing r + I lines zq, «i, . . . , V- These Cl^1^\ minors are called Pliicker coordinates of X. 
They are defined up to a common scalar factor and provide the Pliicker embedding of the grassmannian 
Gl^\ into the projective space P'', where d — C^'^W — 1 (C','J+i[ is the bimomial coefficient). 

The image of Gl^\ in P"* is realized as an intersection of quadrics. This means that the coordinates 
(io, ir) are not independent but obey the Pliicker relations p8| , [ p9| : 

r 

(io,ii, v)(jo, ji, •■•,>) = ^Up^^i^ v)(jo, ■■■jp-i,io,jp+i--,jr) (2.18) 

for all ip,jp, p — 0,1, . . . ,r. Here it is implied that the symbol (zq, ii, . . . , v) is antysymmetric in all 
the indices, i.e., (ip, . . . , ip-i, ip, . . . ,ir) = —(io,---,ip, ip-i, ■ ■ ■ ,ir) and it equals zero if any two indices 



coincide. If one treats these relations as equations rather than identities, then determinants (2.17) would 
give a solution to Hirota's equations. 



The Plucker relations in their general form (2.18) most likely describe fusion rules for transfer matrices 



corresponding to arbitrary Young diagrams. In order to reduce them to the 3-term HBDE, one should 



take ip = jp for p ^ 0,1. Then all terms but the first two in the r.h.s. of (2.18) vanish and one is left 
with the 3-term relation 

{io,ii, . . . , v)(jo, ji,«2, . . . , v) = (jo,«i,«2, • ■ • ,v)(io,ii,i2, • ■ • v) + (ji,«i,i2, • ■ • , v)(jo, «o, «2, ■ • ■ v)- 

(2.19) 
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After the substitution (2.17) these elementary Pliicker relations turn into certain determinant identities. 



For example, choosing = Spj, Xj'J — Sqj, q ^ p, one can recast eq. (2.19) into the form of the Jacobi 



„(i) 



identity: 



D[p\p] ■ D[q\q] - D[p\q] ■ D[q\p] = D[p, q\p, q] ■ D . 



(2.20) 



where D is determinant of a (r+ 1) x (r + l)-matrix and D[pi,p2\qi, 92] denotes determinant of the same 
matrix with pi,2-th rows and (ji,2-th columns removed. Another useful identity contained in eq. (2.1£) 
connects minors D[li, I2] of a (r + 3) x (r + 1) rectangular matrix, where the two rows /i, I2 are removed: 



D[h,h] ■ D[h, h] - D[WM ■ D[l2, h] = D[hM] ■ D[l2,h 



h <h <l3 < h 



(2.21) 



Identifying terms in eq. ( 2. IE ) with terms in HBDE 



one obtains va rious determinant represen- 

(2.22) 



tations of solutions to HBDE. Two of them follow from the Jacobi identity ( 2.20 ): 

ra(/,m) = det (ti(/ + i - a,m - j + a)), i,j = l,...,a, ro(Z,m) = l 

or, in "direct" variables 

T^{u)=det{Tl^,_^{u + i+j-a-l)), z,j = l,...,a, T°(7.) - 1 . (2.23) 

This representation determines an evolution in a from the initial values at a = 1. The size of the 
determinant grows with a. A similar formula exists for the evolution in s: 



T;(w) = det(rf (u + »+ J -,s-l)), i,j = l,...,s, TS{u)^l. 



(2.24) 



The size of this determinant grows with s. Determinant formulas of this type have been known in the 
literature on quantum integrable models (see e.g. [^). They allow one to express Tg{u) through Tf (w) 
or Tl{u). 



A different kind of determinant representation follows from (2.21): 
Ta{l,m) = detAfy, 



0^ 



hi{u + s + a + 2j) if j = 1, k — a; i — 1, k 
hi{u — s + a + 2j) if j = k — a + 1, k; i — 1, k 



(2.25) 



where hi{x) and hi{x) are 2k arbitrary functions of one variable. The size of this determinant is equal 
to k for all < a < fc. This determinant formula plays an essential role in what follows. 

The determinant representations give a solution to discrete nonlinear equations and expose the essence 
of the integrability. Let us note that they are simpler and more convenient than their continuous 
counterparts. 



2.3 Examples of difference and continuous Ai-type equations 

For illustrativ e pu rposes we specialize the Hirota equation to the Ai-case and later study it separately. 
At fc = 2 eq. dO) is 



(2.26) 



Ts{u + l)Ts{u - 1) - T,+i(u)T,_i(u) = ^{u + s)^{u - s - 2) 

with the condition T^i{u) — (here we set Ts{u) = T^{u)). 

This equation is known as a discrete version of the Liouville equation pO|] written in terms of the 
T-function. It can be recast to somewhat more universal form in terms of the discrete Liouville field 



Ts+i{u)Ts-i{u) 



(2.27) 



(u + s)(/)(u - s - 2) 

(see (2.6)), which hides the function <j){u) in the r.h.s. of ( 2.26| ). The equation becomes 

Y,iu - l)Y,{u + 1) = (i;+i(ii) + l)iYs-i{u) + 1) . (2.28) 
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(Let us note that the same functional equation but with different analytic properties of the solutions 
appears in the thermodynamic Bethe ansatz ||3l|.) 

In the continuum limit one should put Ys{u) — 6^^ exp(— (^(a;, t)), u ~ S^^x, s ~ S^^t. An expansion 
in (5 ^ then gives the continuous Liouville equation 

a,V-a> = 2expM. (2.29) 



To stress the specifics of the b.c. (2.15) and for further reference let us compare it with the quasiperi- 



odic b.c. Then the Ai-case corresponds to the discrete sine-Gordon (SG) equation p^ : 

T^+\u) ^ e"A^"r;-i(u - 2), (2.30) 

where a and A are parameters. Substituting this condition into ( |l.7| ), we get: 

T}iu + l)T}iu - 1) - T}^,{u)Tl,iu) = e"A2T°(«)r«(u - 2), (2.31) 

T^iu + 1)T°(^ - 1) - T°+i(u)T°_i(^) = e-'^T}iu)T^iu + 2). (2.32) 
Let us introduce two fields p^'" and (p"''^ on the square {s,u) lattice 

T0(7.) = exp(p^'" + ^^^"), (2.33) 

T}{u + 1) = Ai/2 exp(p^'" - ^'*'"). (2.34) 



and substitute them into (2.31), (2.32). Finally, eliminating p'*'", one gets the discrete SG equation: 

sinh(^'*+i'" + ip'-^'" - ^"'"+1 - ^^'"^1) = Asinh(^"+i^" + ^""i'" + ip''''+^ + ^"^"-^ + a) . (2.35) 

The constant a can be removed by the redefinition lys'*'" — ^ ip^'"^ — ja. 

Another useful form of the discrete SG equation appears in variables (u) ( ^.10 ) . Under condition 
( ^.30 ) one has 

X^+\u) = Xr\u-2), X'Xf+\u+l)X!{u)^l, (2.36) 
so there is only one independent function 

Xl (u) EE xs (u) = -e-"A^i exp ( - 2</j"'" - 2^3"^""^) . (2.37) 

The discrete SG equation becomes ||, |33|, ||: 

(X + xJu + l))(X + Xs(u-l)) 

(1 + \Xs{U + 1))(1 + \Xs{U - 1)) 



In the limit A ^ eq. (2.38) turns into the discrete Liouville equation (2.28|) for Ys(u) — —1 — A ^Xs{u) 



3 Linear problems and Backlund transformations 
3.1 Zero curvature condition 

Consider the square lattice in two light cone variables I and m and a vector function ipa(l,m) on this 
lattice. Let La.a'i^i'ni) and AIa,a'{l,'>n) be two shift operators in directions / and m: 

a' 

Y,Ma^a'(.l,ni)^a'{l,m + l) = Ml,rn). (3.1) 

a" 

The zero curvature condition states that the result of subsequent shifts from an initial point to a fixed 
final point does not depend on the path: 

L{1, m) ■ M{1 + l,Tn) ^ M{l,m) ■ L{l,m + 1). (3.2) 
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HBDE (|Tj) possesses 
matrices: 



a zero-curvature representation by means of the following two-diagonal infinite 



where 



V," 



W = 



= Say-l+Sa.a'Vi" , 
= Sa,a' + Sa,a' + lW^ , 

Ta{l + f,m)Ta+i(^,m) 
Ta{l,m)Ta+l{l -|- l,r7l) ' 

Ta-l{l,m-[- l)Ta+i(/,m) 



(3.3) 



Ta{l,m)Ta{l,m+l) 

More precisely, the compatibility condition of the two linear problems 

1pa{l,m) - ijja+lil + l,m) = Vi^ljjail + , 



(3.4) 



combined with the b.c. (2.14) yields HBDE (1.8). Introducing an unnormalized "wave function" 

fail,rn) = '4ja{l,m)Ta{l,m) , 
we can write the linear problems in the form 

Ta+l{l + l,m)fa{l,m) - Ta+l{l,m)fa{l + l,r7l) = Ta{l , m) f a+l{l + l,m) , 
Ta{l,m + l)fa{l,m) - Ta{l,m)fa{l,m -I- 1) = Ta+lil,m)fa-lil,m+ 1) , 

or in "direct" variables 

T^+l{u)F''{s,u) - T^+\u - l)F''{s + l,u+l)= T^{u)F''+^{s + l,u), 
T^+i{u - l)F''{s, u) - T^{u)F''{s + 1, u - 1) = T^+'^iu - l)F''-\s + l,u), 



(3.5) 
(3.6) 

(3.7) 



(3.8) 



where F°'{1 + mj ~ rn — a) = fa{l, m). 

An advantage of the light cone coordinates is that they are separated in the linear problems (there 
are shifts only of / (m) in the first (second) eq. (Bj)). 



The wave function and potential possess a redundant gauge freedom 

X{a-l) 



X{a - I) 
with an arbitrary function x- 



(3.9) 



The b.c. (2.4) implies a similar condition for the object of the linear problems 

F°-{s, u) = as a < and a> k -1 



(3.10) 



so that the number of functions F is one less than the number of T's. Then from the second equation of 
the pair (^^) at a = and from the first one at a = fc — 1 it follows that F"(s, u) (F'^^^(s, u)) depends 
on one cone variable u + s (resp., u — s). We introduce a special notation for them: 



F"{s,u) = Qk-i{u + s), F''-\s,u)^Qk-i{u-s) 



(3.11) 



Furthermore, it can be shown that the important condition (2.15) relates the functions Q and Q: 

Qk-i{u) = QkMu~k + l). (3.12) 

The special form of the functions F"^ at the ends of the Dynkin graph {a ~ 0, k — 1) reflects the speciflcs 
of the " Liouville-type" boundary conditions. This is to be compared with nonlinear equations with the 
quasiperiodic boundary condition (2.30): in this case all the functions F depend on two variables and 
obey the quasiperiodic b.c. 
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3.2 Continuum limit 



In the continuum limit / = — (5t+, m — ~St^, Ta S'^ Ta, fa ^ <5° ^"/a, S ^ 0, we recover the 
auxihary Unear problems for the 2D Toda lattice {d± = d/dt±): 



Ta 

di Ta-]-lTa—l . 
-Wa = ^ Va-1 , 



(3.13) 



or, in terms of fa, 



Ta+ld+fa - {d+Ta+l)fa = Tafa+l , 
Tad- fa - {d-Ta)fa = Ta+lfa~l 



(3.14) 



The compatibility condition of these equations yields the first non-trivial equation of the 2D Toda lattice 
hierarchy: 

d+Tad-Ta - Tad+d-Ta = Ta+lTa-l- (3.15) 



In terms of 



it has the familiar form 



(Pa{t+,t-) = log 



Ta+l{t+,t-) 
Ta{t+,t-) 



d+d-(pa = e 



(3.16) 



3.3 Duality 

The discrete nonlinear equation has a remarkable duality between "potentials" and "wave functions" 
F"' first noticed in [^ . In the continuum version it is n ot so transparent. Eqs.(B^) are symmetric under 
the interchange of F and T. Then one may treat ( |3.§| ) as linear problems for a nonlinear equation on 
_F's. It is not surprising that one again obtains HBDE (1.7): 



F°(s,u + l)F°(s,u - 1) - F''(s + l,u)F°(s - l,u) = F''+\s,u)F''-\s,u) 



(3.17) 



Moreover, conditions ( 3.10 )-( 3.12 ) mean that even the b.c. for F°'{s, u) are the same as for T^{u) under 
a substitution 4'{u) by Qk-i{u). The only change is a reduction of the Dynkin graph: k ~* k ~ 1. Using 
this property, one can successively reduce the Afc_i-problem up to Ai. Below we use this trick to derive 
Ak-i ("nested") Bethe ansatz equations. 



3.4 Backlund flow 



To elaborate the chain of these transformations, let us introduce a new variable t — 0, 1, . . . , fc to mark 
a level of the flow A^^i Ai and let F^_^-f^{s,u) be a solution to the linear problem at {k — t)-th level. 
In this notation, F^{s,u) = T^{u) and f^_i(s,M) = F '^{s,u) is the corresponding wave function. The 
wave function itself obeys the nonlinear equation (3.17), so F^_2{s, u) denotes its wave function and so 
on. For each level t the function F^{s,u) obeys HBDE of the form (3.17) with the b.c. 



F°(s, u) = as a < and a> t . 



(3.18) 



As a consequence of (3. IS), the first and the last components of the vector F^{s,u) obey the discrete 
Laplace equation (2.3) and under the condition (3.11) are functions of only one of the light-cone variables 
[u + s and u — s respectively). We denote them as follows: 



FO(s, u) = Qtiu + s), Ft\s,u) = Qt{u- s) 



(3.19) 



where it is implied that Qk{u) = 4'{u). It can be shown that ellipticity requirement (ii) and condition 
(2.14) impose the relation Qt{u) = Qt{u — t). 
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In this notation the linear problems (3.8) at level t 



F,%'{s + l,u)Ft{s, u) - ^,7/(5, u - l)Ftis + l,u+l) = F,\,is, u)Fr'is + l,u), 
+ l,u- l)F,-{s,u) - F,\,{s,u)F,-{s + l,u-l)^ F,-^,'(s,u^ l)FrHs + 1,^) 



(3.20) 
(3.21) 



look as bilinear equations for a functions of 4 variables. However, eq. (3.20) (resp., eq. (3.21)) leaves 
the hyperplane u — s + a = const (resp., u + s + a = const) invariant, and actually depends on three 
variable. 



Restricting the variables in eq. ( p.2C| ) to the hyperplane u — s + a = v (where w is a constant), by 
setting 

Tu{t,a) = F^_t{u + a-v,u) (3.22) 
we reduce eq. (3.20) to the form of the same HBDE ( |l.8| ) in cone coordinates t and a. The b.c. is 



Tuit, 0) = Qk-t{'2u - v), T„(t, k-t) ^ Qk-t{v + t-k) = const. 



Similar equations can be obtained from the second linear problem ( 3.21 ) by setting 

f„(6, t) = F^zl'\v + b-u,u + t- k) 
{v is a constant). This function obeys eq. ( |l.8[ ), 

f„ (6+1, t)Tu {b,t + 1) ~ fu{b,t)fu{b + l,t + 1) = fu+i{b + 1, t)fu-i{b, t + I) 
where t now plays the role of the light cone coordinate m. The b.c. is 

Tu(0,t) = Qk-t{2u + t- k-v), fu{k ~t,t) = Qk-t(v) = const. 



(3.23) 

(3.24) 
(3.25) 
(3.26) 



It is convenient to visualize this array of r-functions on a diagram; here is an example for the ^3-case 
(A: = 4): 



1 

Qi{u + s) 

Q2{u + s) 

Qsiu + s) 
4>[u + s) 





Qi{u - s) 
Fiis,u) 
Fi{s,u) 





Q2iu - s) 
Fiis,u) 





Qsiu - s) 
T!iu) 



(3.27) 



{u - s) 



Functions in each horizontal (constant t) slice satisfy HBDE ( |3.17D , whereas functions on the u — s + a ■ 
const slice satisfy HBDE (1.8) with t, a being light cone variables I, m respectively. 



A general solution of the bilinear discrete equation ( ^ ) with the b.c. ( 2.14 ) is determined by 2k 
arbitrary functions of one variable Qt{u) and Qt(u), t = l,...,fc. The additional requirement (ii) of 
ellipticity determines these functions through the Bethe ansatz. 



3.5 Nested Bethe ansatz scheme 

Here we elaborate the nested scheme of solving HBDE based on the chain of successive Backlund trans- 
formations (Sect. 3.4) . This is an alternative (and actually the shortest) way t o obtain nested Bethe 
ansatz equations ( 3.31 ). Recall that the function Tu{t,a) = F^_^{u + a,u) (|^2^) (where we put w = 
for simplicity) obeys HBDE in light cone variables: 



Tu{t + l,a)T„(t, a + 1) - Tu{t,a)Tu{t + l,a + 1) = Tu+i{t + l,a)Tu-i{t,a + 1) 



(3.28) 



Since T„(t, 0) = Qk-t{2u), nested Bethe ansatz equations can be understood as "equations of motions" 
for zeros of Qt{u) in discrete time t (level of the Bethe ansatz). The simplest way to derive them is to 
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consider the auxiliary linear problems for eq. (3.2J 
the simplest possible form. 

Let us assume that Qt{u) has the form 



Here we present an example of this derivation in 



Mt 



Qt{u) 



(3.29) 



(note that we allow the number of roots Mt to depend on t). Since we are interested in dynamics in t 
at a fixed a, it is sufficient to consider only the first linear equation of the pair (3.7): 

T„+i(t + l,a)fu{t,a) - Tu+i{t,a)fu{t + l,a) = r„(t, a)/„+i(i + l,a) . (3.30) 

An elementary way to derive equations of motion for roots of Tu{t,0) is to put u equal to the roots 
of fu{t + 1,0), fu{t,0) and fu+i{t + 1,0), so that only two terms in ( |l30| ) would survive. Combinina; 
relations obtained in this way, one can eliminate /'s and obtain the system of equations 



2)Qt+i(u* 



L(^)gt(^.*+2)Q,+i(u*-2) ^^-^^^ 

as the necessary conditions for solutions of the form ( 3.29| ) to exist. In the more detailed notation they 
look as follows: 

Mt- 



fc=l 



fe=l 



With the "boundary conditions" 



2)) <viu] 



2)) 



Mt+i 

n 



(7(77(1 



,t+i 



)) 



(7(r,(u* - ul + 2)) J-^A a{rj{u] 



,t+i 



-2)) 



Qo(w) 



1, 



(3.32) 



(3.33) 



this system of Mi + M2 + . . . + M^-i equations is equivalent to the nested Bethe ansatz equations 
for j4fe_i-type quantum integrable models with Belavin's elliptic i?-matrix. The sam e eq uations can be 
obtained for the right edge of the diagram ( 3.27 ) from the second linear equation in (^!^). In Sect. 5 we 
explicitly identify our Q's with similar objects known from the Bethe ansatz solution. 



Let us remark that the origin of equations ( 3.32 ) suggests to consider them as equations of motion 
for the elliptic Ruijsenaars-Schneider model in discrete time. Taking the continuum limit in t (provided 
Mt — M does not depend on t), one can check that eqs. ( 3.32 ) do yield the equations of motion for the 
elliptic RS model ||l7| with AI particles. The additional limiting procedure ry — > with finite rjUj = Xj 
yields the well known equations of motion for the elliptic Calogero-Moser system of particles. 

However, integrable systems of particles in discrete time seem to have a richer structure than their 
continuous time counterparts. In particular, the total number of particles in the system may depend on 
(discrete) time. Such a phenomenon is possible in continuous time models only for singular solutions, 
when particles can move to infinity or merge to another within a finite period of time. 



Remarkably, this appears to be the case for the solutions to eqs. (3.32) corresponding to eigenstates 
of the quantum model. It is known that the number of excitations Mt at t-th level of the Bethe ansatz 
solution docs depend on t. In other words, the number of "particles" in the corresponding discrete time 
RS model is not conserved. Though, the numbers Mt may not be arbitrary. 

In the elliptic case degrees of the elliptic polynomials Qt{u) are equal to Mt = N/k (provided 77 is 
incommensurable with the lattice spanned by Wi, bJ2 and N is divisible by k). This fact follows directly 
from Bethe equations (3.31). Indeed, the elliptic polynomial form ( 3.29] ) implies that if u* is a zero of 
Qt{u), i.e., Qt{u*j) = 0, then u* + 2niUJi + 2r7,2W2 for all integers ni,7i2 are its zeros too. Taking into 
account the well known monodromy properties of the cr-function, one concludes that this is possible if 
and only if 

Mt+i+Mt-i^2Mt, (3.34) 

which has a unique solution 

N 

Mt = -t (3.35) 



satisiying b.c. ( 3.33 ). This means that the nested scheme for elliptic Afe-i-type models is consistent 
only if N is divisible by k. In trigonometric and rational cases the restrictions on degrees of QtS are not 
so strong. 



12 



4 The y4i-case: discrete Liouville equation 



In this section we consider the Ai-case separately. Although in this case the general nested scheme is 
missing, the construction is more explicit and contains familiar objects from the Bethe ansatz literature. 



4.1 General solution 

The basic functional relation is 

Ts{u+l)Ts{u 



I) 



Ts+i{u)Ts-i{u) 



where we consider a more general b.c. parametrized by the functions 
auxiliary linear problems (|3.8[) acquire the form 



u + s)(j){u — s), 

1 and set Ts{u) 



T,+i{u)Qiu + s)- Tsiu - l)Q{u + s + 2) = (t>{u + s)Qiu - s - 1), 

T,+i{u)Q{u - s + 1) - Tsiu + l)Q{u -s~l)^4>{u~ s)Q{u + s + 2) 
Here we set Q{u) = Qi{u) and (j){u) ~ Q2{u). Rearranging these equations we obtain 

4){u - 2)Q{u + 2) + <t){u)Q{u - 2) = A{u)Q{u), 

4>{u)Q{u + 3) + 4>{u + 2)Q{u - 1) = A{u)Q{u + 1), 
where we introduce the quantities 

(u - 2)Ts+i{u ~ s) + (j){u)Ts-i{u ~ s - 2) 



Aiu) = 



Aiu) 



Ts{u - s - 1) 
{u + 2)Ts+i{u + s) + 4){u)Ts-i{u + S + 2) 



Ts{u - 



1 



(4.1) 
T}{u). The 

(4.2) 
(4.3) 

(4.4) 
(4.5) 

(4.6) 
(4.7) 



Due to consistency condition ( |4.l| ) A{u) and A{u) are functions of one variable and do not depend on s. 
The symmetry between u and s allows one to construct similar objects which in turn do not depend on 
u. Functions A{u) and A{u), in the r.h.s. of (4.4), (4.5) are the conservation lows of the s-dynamics. 

Running ahead, let us note that the connection between (f) and (f), 4>{u) = (j}{u — 2), and its consequence 
r_i(?i) = (see ( ^.15[ )), simplifies eqs. (4.4)-(4.7). Putting s = and using the b.c. T^i{u) — 0, we find 



A{u) ^A{u) =Ti(w). 

Therefore, the following holds 

Ts{u - l)Ti{u + s) = 4>{u + s - 2)Ts+i{u) + cb{u + s)T,_i(w - 2), 

Ts{u + l)Ti{u - s) = 4>{u - s)Ts+i{u) + 4>{u~s^ 2)T,_i(u + 2), 
(j){u - 2)Q{u + 2) + (f>{u)Q{u - 2) = Ti{u)Q{u) . 



(4.8) 

(4.9) 

(4.10) 
(4.11) 



These equalities are known as fusion relations ||^, |3^ while eq. (4.11) is Baxter's T-Q-relation 

0. So Baxter's Q function and the T-Q-relation naturally appear in the context of the auxiliary linear 
problems for HBDE. 

A general solution of the discrete Liouville equation (for arbitrary and cj)) may be expressed through 
two independent functions Q{u) and Q{u). One may follow the same lines developed for solving the 
continu ous classical Liouville equation (see e.g. [sgj], po| and references therein). Let us consider eq. (4.4) 
(resp., (4.5)) as a second order linear difference equation, where the function A{u) {A{u)) is determ ined 
from th e in itial data. Let R{u) (resp., R{u)) be a second (linearly independent) solution of eq. (4.4) 
(resp., (4.5)) normalized so that the wronskians are 



W{u) 



R{u) Q{u) 
R{u + 2) Q{u + 2) 



(j){u), 



(4.12) 
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R{u + 2) Q{u + 2) 



= 4>iu + i). 



W{u) = 

Then the general solution of the eq. ( |4.1| ) is given in terms of Q and R: 

Ts{u 



Q{u + s + l) R{u + s + l) 
Q{u — s) R{u — s) 



(4.13) 



(4.14) 



This formula is a particular case of the general determinant representation ( p. 251 ) . 

For any given Q{u) and Q{u) the second solution R(u) and R{u) (defined modulo a linear tra nsfor - 



matio n R{u) R{u) + aQ{u) ) can be explicitly found from the first order recurrence relations ( 4.12 ), 
(4.13), if necessary. Let Q{uo) and R{uo) be initial values at u = uq. Then, say, for even r > 0, 



r/2 



(uo + 2j - 2) 



Riuo) 



i?(?io + r) = Q(mo + r) - V — r,^^r,^ \ 

~^ Q("o + 2j)(3(uo + 2j - 2) Q[uq) 



(4.15) 



and so on for other r's and R{u) 



Finally, one can express solution to eq. (4.1) through two independent functions Q{u) and Q{u): 

ro(u-l) , ^ ^{u + 2j-2) 



Tsiu + s - 1) = Q{u + 2s)Q{u - 1) 



Q{u)Q{u 



i=i 



where To(m) can be found from (4.16) by putting s = 0: 

To(w-l) , To{u + l) (j){u) 



Q{u + 2j)Q{u + 2j -2) 



0(u) 



Q{u)Q{u-l) Q{u + 2)Q{u + l) Q{u)Q{u + 2) Q{u - l)Q{u + 1)' 
Note also the following useful representations: 

A{u) = Q{u + 2)R{u - 2) - R{u + 2)Q{u - 2), 
A{u) = R{u + 3)Q(u - 1) - Q{u + i)R{u - 1), 



(4.16) 

(4.17) 

(4.18) 
(4.19) 



which are direct corollaries of (14), (4.5) 



4.2 Equivalent forms of Baxter's equation 



The key ingredient of the construction is Baxter's relation (4.11) and its "chiral" versions (4.4), ([4.5|). 
For completeness, we gather some other useful forms of them. 



Consider first "chiral" linear equations (4.4), (4.5) (thus not implying any specific b.c. in s). Assum- 
ing that Ts{u) obeys HBDE (4.1), one can represent eqs. (4.4), (4.5) in the form 



Ts{u) 




Ts+i{u - 


1) 


Q{u + 


s + 


1) 


Ts+i(u + 


1) 


Ts+2{u) 




Q{u + 


s + 


3) 


Ts+2{u + 


2) 


Ts+3{u + 


1) 


Qiu + 


s + 


5) 


Ts{u) 




Ts+i{u + 


1) 


Q{u- 


5) 




Ts+i{u - 


1) 


Ts+2{u) 






s — 


2) 


Ts+2{U - 


2) 




1) 


Q{u~ 


s — 


4) 



= 0, 



(4.20) 



0, 



(4.21) 



respectively. This representation can be straightforwardly extended to the Afc_i-case (see eqs. (5.37) 
(P8|))- 
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A factorized form of these difference equations is 



^2a„ _ + 2)0(^-1 ) ^^ /^^29„ 0(" + 1) 
0(u)Q(ii + l) 



1) 0. 



(4.22) 



(4.23) 



Here e^" acts as the shift operator, e^^ f{u) = f{u + 1), and X{u) {X{u)) stands for any hnear combi- 
nation oiQ{u), R{u) {Q{u), R{u)). 



Specifying eqs. (4.20), (4.21) to the b.c. T-i{u) — (see (pT^)), we see that both of them turn into 
the equation 



^(-l)Tf (u + a - l)X{u + 2a - 2) = 



(4.24) 



that is Baxter's relation ( 4.1l| ). Furthermore, the difference operator in (4.24) admits a factorization of 
the form (E2|): 



Qju) 
Q{u-2)J ' 



(4.25) 



which is equivalent to the well known formula for Ti(u) in terms of Q{u) 



4.3 Double-Bloch solutions to Bctxter's equation 



In this section we formulate the analytic properties of solutions to Baxter's functional relation (4.4) that 
are relevant to models on finite lattices. 

First let us transform Baxter's relation to a difference equation with elliptic (i.e. double-periodic 
with periods 2LUi/ri, 2uj2/vi) coefficients. 

The formal substitution 

Q{u)P{u) 



with a (as yet not specified ) function P{u) yields 



(t){u - 2) 



ox P{u + 2)<i){u-A) ~ . . A{u)P{u + 2) ^ , . 



P{u~2)<f){u-2) 



(l){u)P{u) 



(4.26) 



(4.27) 



Below we restrict ourselves to the case when the degree N of the elliptic polynomial 4){u) ( |l.3| ) is 
even. Then for any P{u) of the form 



N/2 

P{u) = \{a{i^{u-p,)) 



(4.28) 



with arbitrary pj the coefficients in (4.27) are elliptic functions. Indeed, for the coefficient in front of 
^{u — 2) this is obvious. As for the coefficient in the r.h.s. of (4.27), its double-periodicity follows from 
the "sum rule" (2.. 



Let us represent 4>{u) in the form 

</>(«) = 00 (7/)0i(«), (4.29) 

where (f)oiu), </>!(«) are elliptic polynomials of degree N/2 (of course for N > 2 there are many ways to 
do that). Specifying P{u) as 



Piu) = Mu-2) 



we rewrite (4.27) in the form 



*(u + 2) 



(j)o{u - 4)01 (m) 
0o(u-2)0i(u-2) 



^-(^-2) = 



A{u) 



4>o{u)4>i{u - 2) 



(4.30) 



(4.31) 
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where 



Q{u) 



(4.32) 



Now, the coefficients in eq. (4.31) being double-periodic, it is natural to consider its double-Bloch 
solutions. A meromorphic function f{x) is said to be double-Bloch if it obeys the following monodromy 
properties: 

fix + 2iu^)^B^f{x), a =1,2. (4.33) 

The complex numbers Ba are called Block multipliers. It is easy to see that any double-Bloch function 
can be represented as a linear combination of elementary ones: 



M 



fix) = c,4>(x - X,, z)tel'^, 



i=l 



where [h6l 



^{x, z) 



a{z + X + rj) 



a{z + ri)a{x) 
and complex parameters z and k are related by 



c/{27j) 



Ba = K^'^"/" 



exp(2C(t^Q)(2 -I- 77)) 



a{z + rj) 

a{z - 1]) 



a-{z + T]) 



(4.34) 



(4.35) 



(4.36) 



(C(a;) = g'{x)Iu{x) is the Weierstrass ^-function). Considered as a function of z, <i>(a;, z) is double- 
periodic: 

^{x, z -I- 2u:a) = <&(a;, z). 

For general values of x one can define a single- valued branch of $(a;,z) by cutting the elliptic curve 
between the points z = ±77. In the fundamental domain of the lattice defined by the function 
$(a;, z) has a unique pole at the point 2: = 0: 



^{x, z) 



1 



0(1) 



Coming back to the variable u = x/r;, one can formulate the double-Bloch property of the function 
^(u) ( 4.32| ) in terms of its numerator Qiu). It follows from (4.34) that the general form of Qiu) is 



M 



Qiu) = Qiu; v) = e-'" [] air^iu - u,)) , 
i=i 



(4.37) 



where M = N /2 and v determines Bloch multipliers. 



For the trigonometric and rational degeneration of eqs. ([l.4|), (4.31), (4.37) the meaning of v is quite 
clear: it plays the role of the "boundary phase" for twisted b.c. in the horizontal (auxiliary) direction. 
For each v eq. (4.11) has a solution of the form (4.37). The corresponding value of Ti(u) = Aiu) depends 
on 1/ as a parameter: Ti(m) = Tiiu; v). If there exist v ^ v' such that Ti(m; v) = Ti(u; v'), one may put 
Qiu) = Qiu^v), Riu) = Qiu]v'). In the elliptic case the boundary phase in general is not compatible 
with integrability and so v should have a different physical sense which is still unclear. 



4.4 Bethe equations 

It can be shown that for double-Bloch solutions the relation between (f) and 0, (i>iu) 

Qiu) = Qiu - I), i?(u) = i?(u - 1) , 



so that (see (4.14) 



Tsiu) 



Q(u + s + l) i?(M + S + l) 

Qiu-s-l) Riu-s-l) 



(u — 2), implies 
(4.38) 

(4.39) 
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It is clear that if Q{u) and R{u) are elliptic polynomials of degree N/2 multiplied by an exponential 
function (as in (4.37)), Ts{u) has the desired general form (2.7). 

Under condition ( [4.38[ ) eq. (4.16) yields the familiar result: 



T,{u) = Q{u + s + l)Q{u - s - 1) 



+ 2J-1) 



^ Qiu 



: + 2j + l)Q{u ~ s + 2j - 1) 



(4.40) 



This formula was obtained in 



38| by direct solution of the fusion recurrence relations (O), ( 4.10 ) 



Let Uj and Vj, j = 1, . . . ,M, be zeros of Q{u) and R{u), respectively. Then, evaluating (4.12) at 
u = Uj, u — Uj ^ 2 and u = Vj , u = Vj — 2 we obtain the relations 

Uj)^Q{uj+2)Riuj), (f>iuj-2) = ~Qiuj-2)R{uj), (4.41) 



whence it holds 



iuj - 2) 

(t){Vj) 



R{v, 



[v, - 2) R{vj - 2) 



(4.42) 
(4.43) 



Equations ( 4.42| ) are exactly the standard Bethe equations ( p_.2|). We refer to equations ( 4.43| ) as 
complementary Bethe equations. It is easy to check that eqs. ( 4.42|) ensure cancellation of poles in 
(4.40). A more standard way to derive Bethe equations (4.42), ( 4.43| ) is to substitute zeros of Q{u) 
(or R{u)) directly into Baxter's relation (4.11). However, the wronskian relation (4.12) is somewhat 



more informative: in addition to Bethe equations for Uj, Vj it provides the connection (4.41) between 



them. In the next section we derive the system of nested Bethe ansatz equations starting from a proper 



generalization of eq. (4.12) 



In the elliptic case degrees of the elliptic polynomials Q{u), R{u) (for even N) are equal to M = N/2 
(provided rj is incommensurable with the lattice spanned by Wi, 0^2). This fact follows directly from 



Bethe equations (4.42), (4.43) by the same argument as in Sect. 3.5. 



In trigonometric and rational cases there are no such strong restrictions on degrees M and M of Q 
and R respectively. This is because a part of their zeros may tend to infinity thus reducing the degree. 
Whence M and M can be arbitrary integers not exceeding N. However, they must be complementary 
to each other: M + M = A^. The traditional choice is M < N/2. In particular, the solution Q{u) ~ 1 
(M = 0) corresponds to the simplest reference state ("bare vacuum") of the model. 

We already pointed out that the function Q{u) originally introduced by Baxter (see e.g. ||l| and 
references therein) emerged naturally in the context of the auxiliary linear problems. Let us mention 
that for models with the rational i?-matrix this function can be treated as a limiting value of Tg (u) as 
s ^ CX3 H. Rational degeneration of eqs. (2.7), (4.37) gives 



N 



T,(u)=^n(u-zf ), 



Q{u) ^e^Yliu- Uj) , 



where 



A.. = 



sinh(2i^(s + 1)) 



(4.44) 



(4.45) 



(4.46) 



sinh(2i^) 

(The last expression follows from ( 4.4C| ) by extracting the leading term as w ^ 00.) If the "boundary 
phase" —ii'Tj is real and i' ^ 0, one has from (4.39): 



g(w) = ±2sinh(2zy)e''" lim e 



{2s) 



N-M 



(4.47) 



For each finite s > Ts{u) has N zeros but in the limit some of them tend to infinity. The degenerate 
case v = Q needs special analysis since the limits v ^ and s ^ c» do not commute. 



17 



Another remark on the rational case is in order. Fusion relations ( [4.9[), (Oo|) giv e "Bethe ansatz like" 
equations for zeros of Ts{u) (4.44). Substituting zeros of r^(w ± 1) into (|4.9|),(|4l0|) and using one 
finds: 



sinh(2iy(.s + 2)) Hz'-"'' + s - 1) 



sinh(2!ys) 



is) 



n 



sinh(2t^(5 + 2)) ' 



As) _ (s+l) ^ 
k=l Zj -I- 1 

^ - z[!-'^ + 1 



sinh(2i/s) 



- s - 3) 



n 



(s) (s+l) _ , 

/c=l Zj Zfc i 



(4.48) 
(4.49) 



These equations give the discrete dynamics of zeros in s. They are to be compared with dynamics of 
zeros of rational solutions of cl assic al nonlinear equations , |]l5| . It is an interesting open problem to 
find elliptic analogues of eqs. ( 4.47 )-( 4.49 ). 



5 The Ak-i-case: discrete time 2D Toda lattice 
5.1 General solution 

The family of bilinear equations arising as a result of the Backlund flow (Sect. 3.4), 

F^s, u + \)Fns, u-\)- F^{s + 1, u)F^{s ^l,u)= F^+\s, u)Fr\s, u) 
and the corresponding linear problems, 



(5.1) 

(5.2) 
(5.3) 

(5.4) 

They may be solved simultaneously by using the determinant representation ( 2.25 ). The set of functions 
F^{s, u) entering these equations as illustrated by the following diagram: 



F^^+,\s + 1, u)F,-{s, u) - F;'+,\s, u - l)F^{s + 1, w + 1) = F^+^{s, u)F^+\s + 1, w) , 
F^+,{s +l,u~ l)F,-{s, u) - F,\,{s, u)Ft{s + l,u-l)= F,-+,\s, u - l)Fr\s + 1, u) , 



subject to the b.c. 



F^{s, u) = as a < and a> t 






1 















^2° 





Fi 



(5.5) 



Fi 







F° Fl F^ 

(cf. ( [3.27 )). Functions in each horizontal slice satisfy HBDE (5.1). By level of eq. (5.1) we understand 
the number t. Level is introduced for later convenience. At the moment we do not assume any relations 
between solutions at different levels. 



Determinant formula (2.25) gives the solution to these equations for each level t in terms of t arbi- 
trary holomorphic^ functions h\^\u + s) and t arbitrary antiholomorphic functions h^^^ {u — s). This is 
illustrated by the diagrams: 



1 

















,(3) 



h 



1 



(1) 



;(2) 



(5.6) 



it) 



^Here we call holomorphic (antiholomorphic) a function of u + s (resp.. 
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Then, according to (2.25), the general solution to eq. (5.1) is 

/if+^^(u + s-a + 2) 



At+i) 



{u + s - a + 4) 



;(t+i) 



{u — s + a — t) 



hf^'^\u- s + a-t + 2) 



h[^\u + s-a + 2) 
hi^\u + s-a + A) 



h[^\u + s + a) 
h!i\u - s + a-t) 
h['^\u- s + a-t + 2) 



;(t+i) 



{u- 



a + t) 



a + t) 



, (5.7) 



where < a < i + 1 and the gauge functions Xt i Xt (introduced for normalization ) satisfy the 
following equations: 



x?(u+i)xt("-i)=xr'(«)xr'(^), 
x?(u+i)x?(u-i) = xr'(")xr'(^)- 



(5.8) 



(cf. (2^)). The size of the determinant is t + 1. The first a rows contain functions h'f \ the remaining 



t — a + 1 rows contain h[''\ The arguments of h'f\ h!f^ increase by 2, going down a column. Note 



that the determinant in (5.7) (without the prefactors) is a solution itself. At a = (a = t + 1) it 
is an antiholomorphic (holomorphic) function. The required b.c. ( p.lS| ) can be satisfied by choosing 
appropriate gauge functions x?, X?- 



5.2 Canonical solution 



The general solution (5.7) gives the function T^{u) = F^(s,u) in terms of 2k functions of one variable 



and hl._^. However, we need to represent the solution in terms of another set of 2k functions Qt{u) 



and Qt{u) by virtue of conditions (5.4) in such a way that eqs. ( p.2| ) ( p.3[ ) connecting two adjacent levels 
are fulfilled. We refer to this specification as the canonical solution. 

To find it let us notice that at a = eq. ( p.2| ) consist the holomorphic function Qt{u + s) and a 
function . According to eq. (5.7), is given by the determinant of the matrix with the holomorphic 



entries h!'t\u + .s + 1) in the first row. Other rows contain antiholomorphic functions only, so F^{u, s) 



h[^\u + s + l)'r]i{u — s), where ?7i(u — s) are corresponding minors of the matrix (5/7) at a = 1 



Substituting this into eq. ( p.2\ 
gets relations connecting , 

another boundary a = t + 1. As a result one obtains 

,(1 



and separating holomorphic and antiholomorphic functions one 



h)._^ and Qt{u), Qt+i{u). Similar arguments can be applied to eq. (5.3) at 



hY>{u + s) ^Qtiu + s), h[^\u - s) ^ Qtiu ~ s) 



(5.9) 



and 



Qt+iiu - 2)hfl^{u) = 



2) Qt{u-2) 

(u) Qt{u) 

(u) Qt{u) 

{u + 2) Qt{u + 2) 



(5.10) 
(5.11) 
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where 1 < i < t. Functions x, x iii front of the determinant (pj) are then fixed as follows: 



a > 2, 



(5.12) 



Xtiu) = I n + a - i + 2j - 1) 



a < i - 1 



(5.13) 



It is easy to eheck that they do satisfy eqs. (5^). The recursive relations ( 5.10 ), ( 5.11 ) allow one to 



determine func tio ns fa^ ''' and starting from a given set of Qt{u). These formulas generalize wronskian 



relations ( 4.12| ), ( 4.13 ) to the Afe_i-case. 



Let us also note that this construction resembles the Leznov-Saveliev solution to the continuous 
2DTL with open boundaries. 



5.3 The Bethe ansatz and canonical solution 

The canonical solution of the previous section immediately leads to the nested Bethe ansatz for elliptic 
solutions. 

In this case all functions h\ \ hi are elliptic polynomials multiplied by an exponential function: 



Mi" 



M) 



This implies a number of constraints on their zeros. 

The determinant in (5.10) should be divisible by Qt+i{u — 2) and /i|^-^(u), whence 



(5.14) 



(5.15) 



7(1+1)/ t — l.i\ I t — \A\ 

h\ '{u^ ) _ gt(M^. ') 



(i+l)/ t-\,i 



in'-''' -2) Q*(u*-^'^-2)' 



(5.16) 



(5.17) 



where it*- = m*'^. Furthermore, it is possible to get a closed system of constraints for the roots of Qt{u) 



only. Indeed, choosing u = u*, u = u* + 2 in (5.1C), we get 

Qt+i{u'j - 2)Qt_i(w*) = -Qtiu'j - 2)h?\u]) , 

Qt+iiu])Qt~iiu] + 2) = Qt{u] + 2)hf\u]) . 
Dividing eq. (5. IS) by eq. ( ^.19 ) we obtain the system of nested Bethe equations: 



Qt-i{u]+2)Qt{u]-2)Qt+i{v'^ 
Qt^i{u))Qt{u] + 2)Qt+i(u\-2) 



= -1, 



(5.18) 
(5.19) 

(5.20) 



which coincides with (3.31) from Sect. 3.5. 
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Similar relations hold true for the /i-diagram: 



; (^+1)^-74+1 



-1) 



Qt{ur'^' ~ 1) 



Qt+i{u] + l)Qt-i{u] + 1) = Qt(fl* + 2)hf\v'^ , 
Qt+i{u] - l)Qt-i{u] - 1) = ~Qt{u] - 2)/if )(^) . 
Qt-i{u) + l)Qt{M* - 2)Qt+i(M* + 1) 



Qt-i{u\ ~ \)Qtiu\ + 2)Qt+i{u\ - 1) 



= -1 . 



(5.21) 

(5.22) 

(5.23) 
(5.24) 

(5.25) 



These conditions are sufficient to ensure that the canonical solution for T^{u) (i.e., for i^^(s,u)) 
has the required general form (2.7). To see this, take a generic Q-factor from the product (5.12), 
{ Qt+ i{u — a + 2j))~^ . It follows from ( 5.16| ) that at its poles the j-th and j ' + 1-th rows of the determinant 
( ^.7|) become proportional. The same argument repeated for Q-factors shows that Fj'^j(s, u) has no poles. 

Qt{u — t) leads to condition 

t< s< 



Finally, it is straightforward to see from (5.7) that the constraint Qt{u 
(|t5|) (for 



1 two rows of the determinant become equal) . 
To summarize, the solution goes as follows. First, one should find a solution to Bethe equations 



( 3.31 ) thus getting a set of elliptic polynomials Qt{u), t = 1, . . . , fc — 1, Qq{u) = 1, Qk{u) = 0(u) being 
a given function. To make the chain of equations finite, it is convenient to use the formal convention 
Q-i{u) = Qk+i{u) — 0. Second, one should solve step by step relations ( 5.10 ), ( ^.11 ) and find the 



functions h[^\u), h\''\u). All these relations are of the same type as the wronskian relation (4.12) in the 
Ai-case: each of them is a linear inhomogeneous first order difference equation. 



5.4 Conservation laws 



The solution described in Sects. 5.2 and 5.3 provides compact determinant formulas for eigenvalues of 
quantum transfer matrices. It also provides determinant representations for conservation laws of the 
s-dynamics which generalize eqs. (4.6), (4.7) to the Ak-i-cdse. The generalization comes up in the form 
of eqs. ( 4.20 ), ( 4.21 ) and ( 1.24 ). The conservation laws (i.e., integrals of the s-dynamics) follow from the 
determinant representation (5.7) of the general solution to HBDE. 



Let us consider (C^ -f 1) x (C^ + l)-matrices 



Tls,{s, u) EE Tf+B+B'{n -s + B-B'), B,B' = l, 
TSb' (s, u) - T^^B-B'iu + s + B-B'), B,B' = 1, 



1 



(5.26) 
(5.27) 



where is the binomial coefficient. Let T°[P|i?](s, w) be minors of the matrix (5.26) with row P and 
column i? removed (similarly for (5.27)). 



Theorem 5.1 LetT^(u) be the general solution to HBDE given by eg. (5.7). Then any ratio of the form 

T-[P\R]{s,u) 



ja.R 
^P,P 



,is,u) = 



T^[P'\R]{s,u) 



does not depend on s. These quantities are integrals of the s-dynamics: Ip'p,{s,u) 



minors of the matrix { 5.27 ) give in the same way a complimentary set of conservation law^ 

"^compare with ([4.6|)~4.7). 



(5.28) 



,{u). Similarly, 
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A sketch of proof is as follows. 



Consider the Laplace expansion of the determinant solution (5.7) with respect to the first a (holo- 
morphic) rows: 

r;H = ^Vp(" + s)v5p(^-s) (5.29) 
p=i 

Here P numbers (in an arbitrary order) sets of in dices {pi,p2, . . . ,pa) such that k > pi > p2 > ■ ■ ■ > 
Pa ^ 1, ^Jpiu + s) is minor of the matrix in eq. ( ^.7|) constructed from first a rows and columns pi, . . . ,pa 
(multiplied by xt-ii^ + i^pi^ ~ is the complimentary minor (multiplied by xt-ii''^ ^ ■^))- 

Substitute i?-th column of the matrix ( ^.26| ) by the column vector with components ipp{u + 2B), B — 
1, . . . , + 1. The matrix obtained this way (let us call it (T°-''^'^)b,b') depends on i? = 1, . . . , C£ + 1, 
P = 1, . . . , and a = 1, . . . , k — 1. The "complementary" matrix (T°'^'^) b,b' is d efined by the similar 
substitution of the column vector 'ipp{u + 2_B), i3 = 1, . . . , + 1, into the matrix ( [3.27 ). 



Lemma 5.1 Determinants of all the four matrices introduced above vanish: 

det(r'') = det(f '^) = det(r''^-"-P) = det{f''-'^^^) = 0. 



(5.30) 



functions h^\u + s), h\^' (u — s) in eq. (5.7) 



The proof follows from the Laplace expansion (5.29). ^From this representation it is obvious that 



+ 1 columns of the matrices in ( 5.3C| ) are linearly dependent. This identity is valid for arbitrary 



The conservation laws immediately follow from these identities. Indeed, let us rewrite the determinant 
of the matrix ^'"■•^■^ as a linear combination of entries of the i?-th column: 



det(r'^'^'^) = J2 i-^ f^"iJp{u + 2B')T''[B'\R]{s, u) ^ . 

B' = \ 

Dividing by T°[P'|i?](s, u), we get, using the notation ( 5.28| ): 

^■P'+Vp(" + 2P')- 



(5.31) 



^ {-\)^'rp{u + 2B')r^^;p.{s,u) = 

B' = l,B'^P' 



(5.32) 



The latter identity is a system of linear equations for quantities I'^'p,{s,u), I^'piisju) 



ja,P. 



^P'-i P' -^P'+i P' ■ ■ ■' ^c'^+i P' case of general position wronskians of the func- 



Ta.R 



ja,R. 



tions ipp{u) is nonzero, whence system (5.32) has a unique solution for /pp, (s,u). The coefficients of 

the system do not depend on s. Therefore, Ipp,{s,u) are s-independent too. Similar arguments are 
applied to minors of the matrix (5.27). 



Another form of eq. (5.31) may be obtained by multiplication its l.h.s. by 'tpp{u~2s) and summation 
over P. This yields 

CJ + l 

^ i-lfT^+Biu-s + B)T-[B\R]is,u)^0, (5.33) 

5 = 1 

which is a difference equation for T^{u) as a function of the " holomorphic" variable u + s with fixed 
u — s. 



5.5 Generalized BcLxter's relations 



Equation (5.31) can be considered as a linear difference equation for a function ip°'(u) having l inearly 
independent solutions ipp{u). It provides the Afc_i-generalization of Baxter's relations {i.4), (4.5). This 
generalization comes up in the form of cqs. (4.20), (4.21) and (4.24). 



The simplest cases are a = 1 and a — fc — 1. Then there are k+1 terms in the sum (5.31). Furthermore, 
it is easy to see that 

(5.34) 



Vau) = /i«,(u+l), i^r'(^^K-i(j^ 



22 



Then eq. (5.31) and a similar equation for antiliolomorphic parts read: 

fe+i 



(5.35) 



3 = i 



k+l 



J2i-^Yh^''\^ + 23)r''-'[j\k + u) = 0, (5.36) 
where we put R = k + 1 for simphcity. These formulas may be understood as linear difference equations 



of order k. Indeed, eq. (5.35) can be rewritten as the following equation for a function X{u): 



r,Vi(w + i) T,V2(") 

Tj^,{u + k) T}^,^,{u + k-l) 



Ts+k-ii^ - k + X{u + s + l) 
T,Vfe(i*-fc + 2) X{u + s + 3) 

T}+2k-ii^ + ^) X{u + s + 2k+l) 



0, (5.37) 



This equations h as k s olutions ft-^'lj(u), i ~ l,...,fc. One of them is Qk~i = h)^_'_^{u) (see eq.( ^.9D ) 



Similar equation ( 5.36 ) for the antiholomorphic parts, 



Tt\u) Tallin -I) 

Tallin + I) T^Z,\u) 



Tallin ^k + 2) X{u-s + 2) 



T^I,'{u + k) Ttt-i(^i + fc-l) 



fc-i 



T, 



fc-i 

s-2fe+l 



{u + 1) X{u -~s + 2k) 



= 0, (5.38) 



has k solutions /i^*l]^(w), i = 1, . . . ,k. One of them is Qk-i — h^^lii^)- 



Difference equations (5.37), (5.38) can be rewritten in the factorized form. This fact follows from a 
more general statement. Fix an arbitrary level k and set T^{u) — F^(s,m) (as in Sect. 3). Then for 
each J = 0, 1, . . . , fc — 1 it holds: 



where 



. 7^,_i(u + z-l)I^(u-l)7^(«) 
itj [s, u) 



s-\-i- 



Tltl-2iu I + l)l^(«)I^(u - 1) 



The proof is by induction. At j = eq. (5.39) turns into 



(5.39) 
(5.40) 

(5.41) 

(5.42) 
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This means that F {s,u) does not depend on u + s. Further, 



(see (3.8)). The mductive step is then straightforward. The proof of (5.40) is absolutely identical. 
Now, putting J = fc — 1 we get the following difference equations in one variable: 



(e 



(5.44) 

{s,u + s)){e-^^-+^' ~ R^^:l\s,u + s)) . . .{e"^^-+^^ - Rf-^\s,u + s))Qk-i{^^^ 

(5.45) 

Note that operators e*'^= act only on the coefhcient functions in ( ^.44 ), (5.45). These equations provide 
a version of the discrete Miura transformation of generalized Baxter's operators, which is different from 
the one discussed in the Ref. E^] (see also below). 



Coming back to eq. (5.31) and using relations (5.10), (5.11), one gets after some algebra: 

^l'^{u) = /iJ^V + fc ~ 1) = Qi(u + fc - 1) , (5.46) 

i^l{u)^h^^\u)^Q^{u) (5.47) 
(the proof is straightforward but too lengthy to present it here). 



Then, in complete analogy with eqs. (5.37), (5.38), one obtains from (5.31) the following difference 
equations: 



T^+k-ii^ - fc + 1) X{u + s + k-l) 



r,Vfc(u-fc + 2) X{u + s + k + l) 



T^+2k-ii^ + ^) X{u + s + 3k-l) 



0, (5.48) 



ri_i(u + l) Tl.iu) 



Tl,iu + k) Tl,_,{u + k-l) 



Tlk+i(^-k + l) Xiu-s) 
T^_k{u~k + 2) X{u-s + 2) 

Ti_2,.+i(w + l) X(^.-s + 2fc) 



(5.49) 



to which Qi{u) (resp., Qi{u)) is a solution. The other fc — 1 linearly independent solutions to eq. ( 5.4§| ) 
(resp., ( |3.49[ )) are other algebraic complements of the last (first) line of the matrix in eq. ( ^.7D at a = fc — 1 
(a = 1) multiplied by xtzl (xLi)- 

Further specification follows from imposing constraints ( |3.12 ) which ensure conditions (2.4) forced 
by the usual Bethe ansatz. One can see that under these conditions eqs. ( ^l8|) a n d (|5l9| ) become the 
same. Further, substituting a particular value of s, s = — fc, into, say, eqs. ( ^.48 ), (5.37), one gets the 
following difference equations: 



{-l)''T^{u + a - l)Qi(ii + 2a - 2) , 



(5.50) 



a=0 
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Tf (w - a - 1) 



(u - 2a)(j){u - 2a - 2) 



gfc_i(u-2a) = 



(5.51) 



(we remind the reade r th at '/>(u) = Qk(u)). The latter e quat ion can be obtained direc tly from the 
determinant formula ( p.7| ): notice that under conditions (2^) the determinants in eq. (5/7) become 



minors of the matrix hj^_i{u — 2k + 2j), where i numbers columns running from 1 to fc, j numbers lines 



and runs from to fc skipping the value k — a. Taking care of the prefactors in eq. (5.7) and recalling that 



Qk-iiu), one gets eq. ( ^.51 ). These formulas give a generalization of the Baxter equations 



(y), dU), |ii|). 

At last, we are to identify our Qt's with Qt's from the usual nested Bethe ansatz solution. This is 



achieved by factorization of the difference operator in (5.5C) in terms of Qt(u) 



^(-l)'^Tf(w + a- 1)6^'^° 



a=0 



^2a„ _ Qk{u)Qk-i{u - 2) 
Qk{u -2)Qk-i{u) 



Q2{u)Qi{u-2) 
Q2{u-2)Qi{u) 



Qiiu~2) 



(5.52) 



This formula (proved by a straightforward calculation using inductive argument) coincides with the one 
known in the literature (see e.g. Q). It yields (u) in terms of elhptic polynomials Qt with roots 
constrained by the nested Bethe ansatz equations. They ensure cancellation of poles in the r.h.s. The 
l.h.s. of eq. ( ^.52 ) is known as the generating function for Tf (m); T^(u) for s > 1 can be found with the 
help of determinant formula (2.24). 



6 Regular elliptic solutions of the HBDE and RS system in 
discrete time 

In this section we study the class of elliptic solutions to HBDE for which the number of zeros Mt of 
the T-function does not depend on t. We call them elliptic solutions of the regular type since they have 
a smooth continuum limit. Although it has been argued in the previous section that the situation of 
interest for the Bethe ansatz is quite opposite, we find it useful to briefly discuss this class of solutions. 

It is convenient to slightly change the notation: t''"'(x) = t„(— m, —l),x = urj. HBDE (1.8) acquires 
the form 

T'+i'™(a;)T'^'"+i(a;) - t'+i^"+1(x)t''™(.t) = t^+^'"'{x + r^y^^^+^ix - tj) . (6.1) 
We are interested in solutions that are elliptic polynomials in x, 

M 

{x) = l[a{x-x'.n. (6.2) 



The main goal of this section is to describe this class of solutions in a systematic way and, in particular, 
to prove that all the elliptic solutions of regular type are finite-gap. 



The auxiliary linear problems (3.5) look as follows 



= ^.™(, + + ;;:2?(;p)i|"+;| ^'-'"(x) , (6.3) 



_ ,T,(,m/„A , T''™(a; - 7?)t'+i^'" (a; + 77) , 



^ ' ^ ' T'+i.™(x)T'.™(a;) \ u \ J 

(The notation is correspondingly changed: vl''''"(u?7) = ipu{—m,—l).) The coefficients are elliptic func- 
tions of X. Similarly to the case of the Calogero-Moser model and its spin generalizations ||lj], ||l^ the 



dynamics of their poles is determined by the fact that equations ( |6.3| ), (6.4) have infinite number of 
double-Bloch solutions (Sect. 4). 
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The "gauge transformation" /(x) /(x) — f{x)e°'^ {a is an arbitrary constant) does not change 
poles of any function and transforms a doublc-Bloch function into another double-Bloch function. If Ba 
are Bloch niultiphers for /, then the Bloch muhiphers for / are Bi = Bie^'^'^^, B2 — i?2e^'^'^^, where uji, 
UJ2 are quasiperiods of the cr-function. Two pairs of Bloch multipliers are said to be equivalent if they 
are connected by this relation with some a (or by the equivalent condition that the product B^^i?^'^^ is 
the same for both pairs). 



Consider first eq. (6.3). Since / enters as a parameter, not a variable, we omit it for simplicity of the 
notation (e.g. x^'™ — > x™). 



Theorem 6.1 Eq. (6.0) has an infinite number of linearly independent double-Bloch solutions with 
simple poles at the points x™ and equivalent Bloch multipliers if and only if x™ satisfy the system of 
equations 

r7(xr - x7+^)a(xr ^ x^- ~ r,)a{xT - xj-^ + ry) ^ ^ 
-y a(xr - 2:7+1 - 77)a(x™ - x^ + r;)a(xr - x^f'^) 

All these solutions can he represented in the form 

M 



^-'"(x) = ^ c,(m, z, k)$(x - xr, z)k'I'^ (6.6) 



i=l 



($(x, z) is defined in ( 4-3Sl) ). The set of corresponding pairs {z,k) are parametrized by points of an 
braic curve defined by the equation of the form 



M 



R{n, z) = k'^' + ^ r,{z)K^^-' = . (6.7) 

i=l 

Sketch of proof. We omit the detailed proof since it is almost identical to the proof of the corresponding 
theorem in pq| and only present the part of it which provides the Lax representation for eq. (|6.5D. 



Let us substitute the function 5'™(x) of the form (6.6) into cq. ( |6.3[ ). The cancellation of poles at 
X = X™ — rj and x = x™^^ gives the conditions 



M 



KCi{m, z, k) + \i{m) Cj{m, z, K)<i>(x™ — x™ — 77, z) = , (6-8) 



M 



Ci{m + 1, z, k) = pii{m) Cj{m, z, k)<I>(x™+i — x™, z) , (6-9) 



where 



u-(m)^ ^^^-1"^"' ■Z±2 fell) 



Introducing a vector C(m) with components Ci{m, z, k) we can rewrite these conditions in the form 

iC{m) + K,I)C{m) =0, (6.12) 

C{m + 1) = Mim)C{m) , (6.13) 
where / is the unit matrix. Entries of the matrices C{m) and A4(m) are: 

£y (m) = A,(m)$(xr - xf - 77, z), (6.14) 
M,jim) = /i,(m)$(xr+' - xf, z). (6.15) 
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The compatibility condition of ( 6.12| ) and ( 6.13| ), 

£(m + l)M{m) = M{m)C{m) 



(6.16) 



is the discrete Lax equation. 

By the direct commutation of the matrices Ai (making use of some non-trivial identities for the 
function $(a;,z) which are omitted) it can be shown that for the matrices C and M. defined by eqs. 
(3.14), ( 6.1Cl| ) a nd (|6.15 ), ( |6.11 ) respectively, the discrete Lax equation ( |6.1(^ ) holds if and only if the 

satisfy eqs. ([6.51). It is 
form 



worthwhile to remark that in terms of Xi (m) , /ii (to) equations ( |6.5[ ) take the 
A4to + 1) = -^j(m), i = l,...,M. (6.17) 



Eq. ( |6.12D implies that 



R(k, z) = dei{C{m) + kI) ^ . 



(6.18) 



The coefficients of R{k, z) do not depend on to due to (3.16). This equation defines an algebraic curve 
(3.7) realized as a ramified covering of the elliptic curve. 



Solutions to eq. (6.5) are implicitly given by the equation 

e{Ux^'"' + U+l + U-m + Z) = 0, 



(6.19) 



where the Riemann theta-function Q{X) corresponds to the spectral curve (|6.7[), (6. IS), components of 



the vectors U, U+, C/_ are periods of certain dipole differentials on the curve, Z is an arbitrary vector. 
Elliptic solutions are characterized by the following property: 2uJiU, i — 1, 2, belongs to the lattice of 
periods of holomorphic differentials on the curve. The matrix C{m) — C{1, to) is defined by fixing x''°'"^° , 
io,mo+i ^ i = 1, . . . , M. These Cauchy data uniquely define the curve and the vectors U, U+, U- and Z 



in eq. ( 6.19 ). The curve and vectors U , U+, U- do not depend on the choice of Iq, rriQ. According to eq. 
(3.19), the vector Z depends linearly on this choice and its components are thus angle- type variables. 



The same analysis can be repeated for the second linear problem ( 6^ ) . Now to enters as a parameter 
and we set x''™ x\ for simplicity. The theorem is literally the same, the equations of motion for the 
poles being 



The corresponding discrete Lax equation is 

t{l + l)M{l)^M{l)t{l), 

where| 

tij{l) = K{l)<^{x[ - x] - ■q,z), 

M^,{l)= fi^{l)^x\+^ ^x\-ri,z), 



and 



nM 



A,(0 



A»(0 = 



T^)a{x\ 



v) 



i+i 



r])a{x 



i+i 



Xs'') 

x[ 



V) 



tM 



(6.20) 

(6.21) 

(6.22) 
(6.23) 

(6.24) 
(6.25) 



All these formulas can be obtained from (6.5), ( |6.10| )-(6.15) by the formal substitutions a;f — > x\, 
According to the comment after eq. (|6.19|), the Cauchy data for the Z-fiow x'°'"^°, 



io+i,mo 

versa. 



are uniquely determined by fixing the Cauchy data x' 



lo,mo 



,'o,™o+i 



for the TO-flow and vice 



'A very close version of the discrete L-M pair appeared first in the Ref. as an a priori ansatz 
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7 Conclusion and outlook 



It turned out that classical and quantum integrable models have a deeper connection than the common 
assertion that the former are obtained as a "classical limit" of the latter. In this paper we have tried 
to elaborate perhaps the simplest example of this phenomenon: the fusion rules for quantum transfer 
matrices coincide with Hirota's bilinear difference equation (HBDE). 

We have identified the bilinear fusion relations in Hirota's classical difference equation with particular 
boundary conditions and elliptic solutions of Hirota equation, with eigenvalues of the quantum transfer 
matrix. Eigenvalues of the quantum transfer matrix play the role of the r-function. Positions of zeros of 
the solution are determined by the Bethe ansatz equations. The latter have been derived from entirely 
classical set-up. 

We have shown that nested Bethe ansatz equations can be considered as a natural discrete time 
analogue of the Ruijsenaars-Schneider system of particles. The discrete time t runs over vertices of the 
Dynkin graph of Ak-i-type and numbers levels of the nested Bethe ansatz. The continuum limit in t 
gives the continuous time RS system This is our motivation to search for classical integrability 

properties of the nested Bethe ansatz equations. 

In addition we constructed the general solution of the Hirota equation with a certain boundary 
conditions and obtained new determinant representations for eigenvalues of the quantum transfer matrix. 
The approach suggested in Sect. 5 resembles the Leznov-Saveliev solution to the 2D Toda lattice 
with open boundaries. It can be considered as an integrable discretization of the classical VF-geometry 

il- 

We hope that this work gives enough evidence to support the assertion that all spectral characteristics 
of quantum integrable systems on finite ID lattices can be obtained starting from a classical discrete 
soliton equations, not implying a quantization. The Bethe ansatz technique, which has been thought 
of as a specific tool of quantum integrability is shown to exist in classical discrete nonlinear integrable 
equations. The main new lesson is that solving classical discrete soliton equations one recovers a lot of 
information about a quantum integrable system. 

Soliton equations usually have a huge number of solutions with very different properties. To extract 
the information about a quantum model, one should restrict the class of solutions by imposing certain 
boundary and analytic conditions. In particular, elliptic solutions to HBDE give spectral properties of 
quantum models with elliptic i?-matrices. 

The difference bilinear equation of the same form, though with different analytical requirements, has 
appeared in quantum integrable systems in another context. Spin-spin correlation functions of the Ising 
model obey a bilinear difference equation that can be recast into the form of HBDE [Q, [^ . 
More recently, nonlinear equations for correlation functions have been derived for a more general class 
of quantum integrable models, by virtue of the new approach of Ref. |^9[| . 

Thermodynamic Bethe ansatz equations written in the form of functional relations |30|, js^ (see 
e.g., 1^^) appeared to be identical to HBDE with different analytic properties. 

All these suggest that HBDE may play the role of a master equation for both classical and quantum 
integrable systems simultaneously, such that the "equivalence" between quantum systems and discrete 
classical dynamics might be extended beyond the spectral properties discussed in this paper. In particu- 
lar, it will be very interesting to identify the quantum group structures and matrix elements of quantum 
L-Operators and i?-matrices with objects of classical hierarchies. We do not doubt that such relation 
exists. 
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